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Abstract 

The microscopic spectral correlators of the Dirac operator in three-dimensional 
Yang- Mills theory coupled to fundamental fermions and with three or more colours 
^ • are derived from the supersymmetric formulation of partially quenched effective 

Lagrangians. The flavour supersymmetry breaking patterns are appropriately iden- 
tified and used to calculate the corresponding finite volume partition functions from 
Itzykson-Zuber type integrals over supersymmetric cosets. New and simple deter- 
minant expressions for the spectral correlators in the mesoscopic scaling region are 
thereby found. The microscopic spectrum derived from the effective finite volume 
partition function of three-dimensional QCD agrees with earlier results based on 
the unitary ensemble of random matrix theory and extends the corresponding cal- 
culations for QCD in four dimensions. 



1 Introduction and Summary 

Quantum chroinodynamics in three spacetinie dimensions (QCD3) provides an interesting 
and sometimes solvable testing ground for phenomena which occur in its four dimensional 
counterpart. It is related to the high temperature limit of QCD4 [^ and also to quantum 
antiferromagnetism [Q. In this paper we will study some aspects of the spontaneous 
breaking of flavour symmetry in three dimensional QCD ||ll,[@]-0- This mechanism is 
the analog of the breakdown of chiral symmetry in four dimensions which is believed to 
be an important property of strong interactions. 

QCD3 with massless quarks in the fundamental representation of the SU{Nc) gauge 
group, with Nc > 3, possesses a continuous, global flavour symmetry group U{Nf) which 
acts on Nf species of two-component, parity-odd complex fermion fields tpi as t/'j 1— > f/j ■^ ipj, 
ip \-^ Uj ip , where U E U{Nf). Symmetry breaking occurs when there is an even 
number Nf = 2nf of fermion flavours . To understand the mechanism, one introduces 
a small fermion mass term iip M.ip with mass matrix M. which preserves the parity 
symmetry of the massless Euclidean field theory but which explicitly breaks the flavour 
symmetry. This can be achieved, for example, by arranging the quark masses mi into 
pairs of opposite sign in the diagonal mass matrixp 

A^ = diag(mi,...,mn^,-mi,...,-m„J . (1.1) 

In this case, the corresponding fermion determinant is positive definite and one may invoke 
the Vafa-Witten theorem to argue that, if flavour symmetry is spontaneously broken, 
then the diagonal elements of the Nf x Nf Hermitian fermion condensate matrix 



E^. = (0 



V'Vj 



(1.21 



are equal in magnitude and of the same signs as the corresponding masses. This implies 
that the global flavour symmetry group is broken according to 



U(2nf) -^ U{nf) x U{nf) (1.3) 



by the order parameter 



2nf 



tr 



;i-4) 



In this instance, the discrete symmetry group Z2, generated by the product of three- 
dimensional spacetime parity and the flavour exchange ipi <-^ 'ipn.f+i, i = l,---,"^/, is 
unbroken and remains a good symmetry even at the quantum level. 



^Recall that in three Euchdean spacetime dimensions the (Dirac) fermion mass term is purely imagi- 
nary and odd under parity, ^ij) 1— > —tpip. Fermion masses may therefore be positive or negative in three 
dimensions and change sign under parity. 



On the other hand, by an analog of the Banks-Casher relation 0, the distribution of 
small eigenvalues of the three-dimensional Euclidean Dirac operator ilp is related to the 
condensate ( |1.4| ) by 

So = -^^ , (1-5) 

where p(A; mi, . . . ,mAr^) is the spectral density of iJj) and V is the volume of three- 
dimensional spacetime. Understanding the function p is therefore tantamount to a de- 
tailed description of the dynamics underlying flavour symmetry breaking in QCD3. This 
distribution is difficult to compute in general, even near the spectral origin A = 0. How- 
ever, since only the low momentum modes of the Dirac operator spectrum are relevant, 
one may propose that p could be computed in the ergodic regime where the zero momen- 
tum mode of the corresponding Goldstone field f/ dominates [^. In this case, the effective, 
finite-volume partition function becomes remarkably simple, and it is equivalent to the 
representation of QCD3 in terms of microscopic degrees of freedom precisely in the limit 
of zero momentum. The spacetime integration over the effective Lagrangian produces an 
overall volume factor V ^ and the partition function simplifies to a finite dimensional group 
integral of the zero modes of V over the coset space determined by the pattern of fiavour 
symmetry breaking. For the symmetry breakdown ( |1.3|) , we have 

where T^ = 1„ ® cts (!«/ denotes the nf x Uf identity matrix and a^ the 2x2 diagonal 
Pauli spin matrix), DU denotes the invariant Haar measure on the Nf x Nf unitary group 
U{Nf), and 

is the corresponding Goldstone manifold. The beauty of the expression ( |1.6| ) comes from 
the observation |^ that the integration may be extended from the symmetric space ( |1.7| ) 
to the group manifold of U{2nf). This follows from the fact that the subgroup of the 
flavour symmetry group whose adjoint action leaves the matrix Fs invariant is precisely 
Uijif) X Uijif), so that the two integrals agree up to the volume of this stability subgroup 
in the Haar measure DU . The finite volume partition function (|1.6| ) may then be evaluated 
analytically using the Itzykson-Zuber formula [0, and thereby used to explicitly derive 
quantities such as spectral sum rules for QCD3. 

One is ultimately interested in taking the thermodynamic limit l^ -^ 00 and the quark 
masses rrii -^ 0. In the ergodic regime, one must keep the linear dimension V^'^ of the 
system much smaller than the Compton wavelength of the Goldstone bosons, which is 
tantamount to holding fixed the parameters [Q 

Ui = VT.Qmi . (1.8) 



This approximation ensures that the non-zero momentum modes factorize from the effec- 
tive Euchdean QCD3 partition function. The crucial observation made some time ago |^ 
(see |Tl| for a recent review) was that the effective partition function ( |1.6| ) can be equiva- 



lently described by the large A^ limit of an A^ x A^ unitary random matrix ensemble. The 
simplest matrix model of this type is defined by the partition functiona 

2 Nf 

DT e-^''^' Y[det{T-tmj) (1.9) 

of the Gaussian unitary ensemble, where DT is the Gaussian-normalized Haar measure 
on the Lie algebra u{N) oi N x N Hermitian matrices. Since the massless Euclidean 
Dirac operator ip in an arbitrary background field is Hermitian, the matrix model ( |1.9| ) 
possesses the same global symmetries as the original field theory. The spacetime volume 
V translates directly into the size A^ of the matrices. The main conjecture put forward 
in this context is that the Dirac operator spectrum can be computed from the spectral 
correlation functions of the matrix model (|1.9| ). With the assumption that the spectral 
properties of random matrix theory carry over to QCD3, the quantity tt/HqV = l/p(0) 
is the mean level spacing between the smallest eigenvalues of the Dirac operator, where 
p(A; rrii, . . . , rriisfj:) may now be computed as the spectral density of the Hermitian matrix 
T in the ensemble (|1.9|) . The ergodic limit of QCD3 thus becomes the microscopic or 
local scaling limit of the random matrix model ( |1.9| ), i.e. N ^ 00 with Nrrii fixed. This 
special limit motivates the introduction of a microscopic spectral density defined in the 
mesoscopic scaling region by 

p.(«;.„...,c.^,) = Jhn^^p(^;^,...,^) , (1.10) 

where u = 7rp(0)A are the unfolded Dirac operator eigenvalues. For broken flavour sym- 
metry the quantity ( p.. 10 ) is a non-trivial function, and the order parameter (|1.4| ) may be 



computed from the resolvent function of the ensemble ( p..9|) as 

1 d 






for any i = 1, . . . ,nf. 

The microscopic spectral density ( |1.10|) has been computed using random matrix the- 



ory techniques in |^ |T^ and related aspects of QCD3 within this framework are described 
in |T3|-|T^. However, in order to have a better understanding of the relationship between 
the effective field theory ( |1.6| ) and the random matrix theory ( |1.9| ) in the microscopic 
domain, one would like to compute the spectral density directly from the low-energy ef- 
fective field theory. Indeed, the calculation of the Dirac operator spectrum should not 



^The universality of random matrix theory results, i.e. the insensitivity to the details of the particular 
matrix potential in the appropriate limit, is well established |12[ . In this paper we will consider the 
simplest Gaussian potentials, consistent with only the general symmetries of the problem as input. 



rely solely on random matrix theory techniques, and should follow directly from quan- 
tum field theory. This is a non-trivial computation because the usual infrared limit of the 
QCD3 partition function, which is dominated by the Goldstone modes associated with the 
spontaneous flavour symmetry breaking, does not access the Dirac operator spectrum. If 
possible though, the matching of such results with those of random matrix theory would 
constitute a direct proof that the universal matrix model calculations do indeed reproduce 
the microscopic Dirac operator spectrum of QCD3. 

The problem was solved for four- dimensional QCD in ||l^ by the introduction of a 



species of fictitious "valence" quarks, paired with yet another set of fictitious particles of 
opposite quantum mechanical statistics. This leads to a model of "partially quenched" 



QCD3 [jT8[ containing A^^ valence quarks and their supersymmetric partners, of masses /Xj 
and /2j, i = 1, . . . , iV„, respectively, and Nf unquenched (physical) sea quarks of masses 
rrii, i = 1, . . . , Nf. The quark fields are all assumed to transform in the fundamental 
representation of the gauge group. In the original field theory formulation, the Euclidean 
partition function is 

Zn„nX{^&, {/i.,/^.}) = JldA] n "^^fZZT-l ft det(.^ - ^m,) e'^-IA] ^^^^^^ 



i=i 



where S'ym[^] is the three-dimensional Yang-Mills action. When Hi = fli for each i = 
1, . . . ,Ny, the fermion determinants arising from integration over the valence quarks are 
cancelled by the contributions from the corresponding bosonic ghost quarks of the same 
masses. Then, the partition function (|1.12| ) reduces to that of ordinary QCD3 with Nf 
physical flavours of fermions. We shall refer to this case as the "supersymmetric limit" , 
but we will keep the masses generically distinct to lift the degeneracy between the valence 
quarks and their superpartners. The partition function ( |1.12| ) is now also the generating 



function for mass- dependent condensates of the extra quark species as |]T9| 

id f ( ^ 



^s{ifj'j;^i,---,ujNf) = ~j;^g—,^^z^f<NA\j;^r,{f^i,f^i} 



{^i=iii} 



= -2^So/.,yrfn^^^-^^^,^^^, (1.13) 

where the second equality follows from the spectral representation of the condensate 
and the fact that ps is an even function of u. In this equation ps is the microscopic 
spectral density of the Dirac operator in the original, unquenched field theory. The 
condensate in ( |1.13|) can be expressed as the Stieltjes transform of ps{u;u!i, . . . ,ujj^j,) 



which, under suitable convergence criteria for the function S^, has a unique inverse given 
by the discontinuity of ( |1.13|) across the real axis. 



^0 






;i.i4) 



Therefore, a detailed understanding of the partially quenched partition functions ( |1.12|) 
will enable a precise, field theoretical determination of the microscopic spectral density 
of the Dirac operator for QCD3. 

In this paper we will discuss how to evaluate the microscopic spectrum of the QCD3 



Dirac operator using the partially quenched quantum field theory ( |1.12| ), thereby extend- 
ing the computations in four dimensions |1^ . In doing so, we will uncover some subtleties 
concerning the breaking of flavour symmetry in these models. The main observation we 
shall make may be summarized as follows. The basic flavour symmetry of the QCD3 
action with Nj sea quarks and A'"^ valence quarks is parametrized by the Lie supergroup 
GL{Nf + Ny\N^). As we are ultimately interested in studying the microscopic density of 
Dirac operator eigenvalues corresponding to a broken symmetry phase, we assume that 
Nf = 2nf and work with the parity-symmetric mass matrix ( |1.1| ). Let us further assume 
that there are n^ positive masses /ij and n~ negative ones, so that N^ = n^ + n~ . We will 
see that in this case the flavour supersymmetry in the massless limit is broken according 
to 



GL{2nf + nt + n; 



\nt + n: 



GL{nf + n^\n^) x GL{nf + 



njn. 



;li5) 



The details of the symmetry breaking pattern ( 1.15 ) present some subtleties in the 
computation of the spectral density via the expression ( [1.13| , prT^ . The low momentum, fi- 
nite volume partition function corresponding to the field theory ( |1.12| ) is a supersymmetric 
generalization of the Itzykson-Zuber integral ( |l.(j| ) taken over the Goldstone supermanifold 



G{nf]n. 



V ' ^v 



GL{2nf + nt + n. 



V \'^v ~^ "^v 



GL{nf + n+\n+ 



X GL{nf + n^ 



;i.i6) 



However, in contrast to the unquenched case, one cannot extend the coset ( |1.16| ) to the full 
flavour supergroup in a straightforward way, because the volume of the unitary supergroup 



vanishes in its Haar-Berezin measure |2T[. Therefore, the finite volume partition function 
in this case must be dealt with as an integral over a coset superspace, rather than a 
Lie supergroup. This makes the parametrization of the integration variables far more 
intricate, and there is no known generalization of the Itzykson-Zuber formula for such 
superspaces (The generalization of the Itzykson-Zuber formula for the unitary supergroup 



has been derived in |^). These subtleties do not arise in the four dimensional case, as 
the patterns of chiral symmetry breaking in both quenched and unquenched cases are 
completely analogous, and the effective finite volume field theory is given by an integral 
over a Lie supergroup which is the straightforward supersymmetric generalization of that 



for the unquenched case fl^ . 

For example, consider the case of only a single species of valence quarks, A^„ = 1, which 
is the pertinent partially quenched model from which to extract the spectral density. Then, 
according to ( p..l5| ), the flavour symmetry breaking pattern is 



GL{2nf + 1|1) — > U{nf) x GL{nf + 1|1) 



;i.l7) 



It follows that there is no symmetry breaking associated with the flavour supersymme- 
try of the theory, only that which is associated with the original unquenched field theory. 
While the symmetry breaking pattern ( |1.17[ ) complicates the evaluation of the function ps 
from (|1.13|jl.l^), it will turn out to be the correct answer which gives the spectral distri- 



butions in the microscopic scaling limit that are anticipated from random matrix theory. 
This will provide an analytical demonstration that the microscopic distribution of eigen- 
values of the QCD Dirac operator in three dimensions can be computed from an intricate, 
supersymmetric extension of the effective finite- volume QCD3 partition function ( |1.6|) , 
i.e. that the Dirac operator spectrum in three dimensions can be derived directly from 
quantum field theory. In addition, it yields an analytic proof that the smallest eigen- 
values of the QCD3 Dirac operator are correlated according to a random matrix model 
whose form is dictated by the global symmetries of ip. In this way we will present the 
appropriate generalization of the results for four spacetime dimensions and the chiral uni- 
tary ensemble of random matrix theory to three spacetime dimensions and the ordinary 
unitary ensemble. 

The organization of the remainder of this paper is as follows. In section 2 we present 
some field theoretical arguments for the symmetry breaking patterns (|1.15|) . In section 3 



these same patterns are derived using a random matrix theory representation of the quan- 
tum field theory (|1.12|) , along with the finite volume, low energy effective field theory in 



the local scaling limit. In section 4 we illustrate these formal properties of the partition 
functions ( |1.12| ) by performing some explicit calculations in the quenched approximation. 
In section 5 we present the calculation of the microscopic spectrum of the Dirac operator. 
There we derive new expressions for the spectral density ps which agree with those pre- 
viously derived in the literature using random matrix theory, but which are much more 
compact and useful. As an interesting by-product of this analysis, we will also uncover 
an elegant representation of the finite volume partition function (|1.6|) itself. In section 
6 we extend this analysis to compute all microscopic /c-point correlation functions. Two 
appendices at the end of the paper contain some technical aspects of our analysis. In 
appendix A we formally prove that the low-energy effective field theory reduces exactly 
to (|1.6|) in the supersymmetric limit, as it should. In appendix B we present a simple and 



self-contained derivation of the Itzykson-Zuber formula for the unitary supergroup [^ 
which is used in the calculations of sections 5 and 6. 



2 Flavour Symmetry Breaking in Three Dimensions 

In this section we will begin our analysis of the microscopic regime of QCD3 by presenting 
some heuristic, field theoretical arguments for the patterns of flavour symmetry breaking 
in the partially quenched models ( |1.12| ). In this paper we will deal only with the case of 



an even number Nf of physical fermion flavours. For an odd number of physical flavours, 
the discrete Z2 symmetry of the theory, composed of parity and flavour exchange, is 



broken explicitly in the massive case, while for massless quarks it is broken radiatively 
by a gauge invariant anomaly which manifests itself in the appearence of a Chern-Simons 
term at one- loop order [Q. For even Nf this anomaly vanishes, and so we shall henceforth 
work with this case to facilitate some of the arguments which follow. 

Consider partially quenched, Euclidean QCD3 with Nc > 3 colours, Nf flavours of 
fundamental sea quarks, and N^ flavours of fundamental valence quarks. As mentioned in 
the previous section, the tree-level global symmetry group of this quantum field theory is 
GL{Nf + Ny\N^). This group rotates the fermionic sea and valence quarks, and also the 
bosonic superpartners of the valence quarks, among each other. To define the quantum 
path integral, we need to choose a measurable subspace of it. The maximally symmetric 
Riemannian submanifold of GL{Nf + A^„|A^^) is supported by the ordinary, compact Lie 
group U{Nf + A^^,) in the fermion-fermion sector, and by the non-compact Lie group 
GL{Ny, C)/U{Ny) in the boson-boson sector. While the Grassmann integrations are well- 
defined in the path integral representations ( |1.12| ), convergence of the integrations over the 



bosonic quark fields is inconsistent with compact U{Ny) flavour rotations in this sector. 
For this reason, the bosonic valence quarks must transform under a non-compact group 
of flavour transformations which is consistent with convergence requirements. We shall 
meet this requirement again in a somewhat more explicit form in the next section. This 
structure is necessary to produce a positive definite quadratic form for the kinetic and 
mass terms of the low-energy effective Lagrangian. 

We will begin by deducing the symmetry breaking pattern (|1.15|) for the massless 
quantum field theory by employing a generalization of the Coleman- Witten argument 
for ordinary QCD4 p3[. For this, it is instructive to first deduce the pattern (pTSp in 
the original, unquenched massless field theory, a possibility which was mentioned in |p. 
Under a global rotation in flavour space, the Hermitian fermion condensate matrix ( |L2| ), 
which is a natural order parameter for flavour symmetry breaking, transforms as 

Si — yU^EU , UeU{Nf), (2.1) 

while under a Z2 parity transformation it maps as 

Si — > -S . (2.2) 



We assume that all of the criteria of ||23[ apply in our case. In particular, the effective 
potential Vj, obtained by integrating out the Yang-Mills fields, is a U{Nf) x Z2 invariant 
function of the condensate matrix S which does not display any accidental degeneracy 
with respect to the flavour and parity symmetries of the theory. This means that any 
ground state of the quantum field theory, obtained by minimizing the effective potential 
V/(S), can be obtained from any other one by a U{Nf) x Z2 transformation. 



From the continuous symmetry (2J^) it follows that Vj is a function only of the Nj real- 
valued eigenvalues o"i, . . . , aNf of the condensate matrix S. From the discrete symmetry 



( p.2| ), it follows that the effective potential is an even function of the eigenvalues, 

Vf (-0-1, . . . , -aNf) = V/(o-i, . . . , aNf) . (2.3) 
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The simplification we would now like to make is to take the limit Nc ^>- oo oi a large 
number of colours. It is well-known that in this limit only the contributions from planar 
't Hooft diagrams survive |]2^. These graphs correspond to the Feynman diagrams that 
contain only a single quark loop which, in the expansion of the invariant function Vf 
in powers of traces of powers of the condensate matrix, are generated by single trace 
insertions of powers of S, i.e. 

V>(S)= trV;(S) + 0(i-)=^V;(a,) + 0(^) , (2.4) 

where V/ is some scalar function which is independent of Nc. Since the eigenvalues ai are 
independent variables, the ground states are determined by minimizing each term V/(o"j) 



in (2^). Thus each eigenvalue of S is at the minimum of the function V/. From the 



reflection symmetry ( p.3| ) it follows that if ctj is a minimum, then so is — ctj. We conclude 



that either all eigenvalues vanish, or else there are Uf = Nf/2 strictly positive and equal 
eigenvalues ai with the nj other ones being their Z2 reflections. The first possibility, 
which corresponds to the case of unbroken flavour symmetry, may be excluded by arguing 
similarly to [^ . Namely, the three-current Green's function contains only massless poles 
in the large Nc limit and so the fermion bilinear current must create a massless scalar 
particle from the vacuum state. By Goldstone's theorem, this implies the spontaneous 
breakdown of the continuous symmetry, and so only the second possibility remains. In 



this way we deduce the symmetry breaking pattern ( |1.3|) . 

Now let us consider the partially quenched quantum field theory ( |1.12|) with massless 
fields. We define a condensate matrix E similarly to (|1.2|) . It is a supermatrix which lives 
in the Lie superalgebra of the flavour supergroup GL{Nf + Ny\Ny) and which transforms 
under the adjoint action of global flavour superrotations. Since the real-valued eigenvalues 
of S are commuting variables, we can go through the above argument by considering 
individually the fermion-fermion and boson-boson blocks of the theory. The total, effective 
potential for the partially quenched model may then be written as 

Nf iv„ N^ _ _ . 1 . 

v(p) = J2vM) + T.^.i^^) + T.Ma) + o(—) , (2.5) 

where V/ is, as above, the contribution from Feynman diagrams consisting of only a single 
physical quark loop, V^ is the contribution from single valence quark loops, and V^ comes 
from single loops of the bosonic ghost quark fields. The ground state is now determined 
by independently minimizing each of the terms Vf{ai), Vv{C,i) and V^i^a) in ( |2.5| ). The 
first term is minimized as explained above. 

For the last two terms, we recall that for our purposes the valence quarks are fictitious 
particles and we are really interested in the near massless limit of the quantum field 
theory (|1.12|) . Let us assume that there are n^ positive masses fii and n~ negative 
ones. Since we are ultimately interested in taking the supersymmetric limit /ij = /ij, 
the same structure is assumed to hold for the fermionic valence quark masses /Xj. Again 

9 



convergence requirements restrict the flavour symmetry of the ghost quark flelds to a 
non-compact U{n^,n~) subgroup of GL{Ny,C), which deflnes a maximally symmetric 
Riemannian submanifold. Let n„ = mm{n^,n~). Since the bosonic ghost quark flelds 
transform as scalars under spacetime parity, it follows that the action in ( |1.12D has a 
"reduced" Z2 parity symmetry deflned by reflecting ny of the masses of a given sign into 
n^ masses of the opposite sign and interchanging the n„ pairs of corresponding flavours. 
In terms of the eigenvalues of S this is represented as the symmetry 



eo 



< 



«+n„ 



1 < a < Uy 



(2.6) 



of the effective potential (pl5| ) (We have used the fact that the eigenvalues in the large Nc 
limit may be arbitrarily ordered). Applying the above reasoning to the minima of Vy we 
conclude that either all eigenvalues C,a vanish in the ground state, or else there are Ny — 2ny 
non-vanishing equal eigenvalues, ny of which map onto the remaining Uy eigenvalues by 
a Z2 parity transformation. Excluding the flrst possibility, this breaks the Uiriy^riy) 
flavour symmetry to the subgroup U{ny) x U{n~). By super symmetry, we may also infer 
the symmetry breaking U{Ny) — > U{Ny — Uy) x U{ny) for the valence quarks. Since 
the pattern of symmetry breaking is determined entirely by the pattern of eigenvalues at 
the minimum of the effective potential, we arrive at the flavour supersymmetry breaking 
pattern (|1.15|) . 

Generally, by a supersymmetric generalization of the Vafa-Witten theorem [^] and the 
appropriate assumption of spontaneous flavour supersymmetry breaking, one may arrive 



at the pattern ( 1.15 ). The low energy physics of QCD3 is completely determined by the 
spontaneous symmetry breaking pattern described above p[ . The effective fleld theory for 
the low momentum modes of the Goldstone superfleld may now be derived analogously to 
the four dimensional case which utilizes chiral perturbation theory [IT, |TB[. In the ergodic 
regime, the effective Lagrangian in Euclidean space is given by 



£ 



eff 



^ STi OMdM-^ STr 

4 '^ ^ 2 



M' 



u + u- 



(2.7) 



where /,r is the pion decay constant, AA' is the mass matrix of the fleld theory ( |1.12| ), and 
the supertrace STr will be deflned in the next section. The masses of the Goldstone bosons 
are given by the usual Gell Mann-Oakes-Renner relation M\q = {M-aa + ■M.bb)'^o/ f^- 
The corresponding superflelds U live in the vacuum manifold for the symmetry breakdown 
and may be parametrized as 



W = [/diag(l„^ 



+n+; ^r 



.f+r. 



1. 



-1. 



u 



-1 



(2. 



with U G Q{nf;ny,ny). This leads to a supersymmetric generalization of the flnite 
volume partition function ( |1.6D with integration domain the symmetric superspace ( |1.16|) . 

A tantalizing aspect of the symmetry breaking here is the possibility of unbroken 
flavour supersymmetry in the valence sector, which occurs whenever Uy = 0, i.e. all valence 
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quark masses are either positive or negative, thereby destroying the parity symmetry of the 
valence fields. Such a situation arises, for example, in the fully quenched case Nf = with 
a single species of valence quarks. Then, the GL{1\1) fiavour supersymmetry is unbroken. 
In that case, a very simple argument ^ is sufficient to determine the precise form of the 
partition function ( |1.12| ) in the microscopic limit. For this, we expand the vacuum energy 
in powers of the masses fi and /i by treating the mass terms as perturbations. Since there 
is no spontaneous breaking of the continuous fiavour symmetry, there are no massless 
particles in the spectrum of the quantum field theory, and hence the perturbation series 
does not produce any infrared divergences. The vacuum energy per unit spacetime volume 
thereby admits the Taylor series expansion 

-—lnZo^l{H,fl) = Zo + ZifX + Zifj,+ Yl Znm fJ'"' P'"^ ■ (2.9) 

■'^ n+m>2 

The constant zq affects only the overall normalization of the fully quenched partition 
function and may be set to zero. In the microscopic limit, the masses are taken to vary 
as yU, /i ~ -^. In the thermodynamic limit, the infinite sum in ( |2.9| ) therefore vanishes and 
we have 

-lnZo,i{fi,fi) = Nz,fi + Nz^^ + (j^^ . (2.10) 



The constants Zi and Zi may be determined from the definitions ( |1.4| ) and ( |1.13| ) to be 

1 9 In Zo,i (/i,/i) 



Zl 



N dfi 

1 glnZo_i(/i,/i) 



= +zSo sgn/i , 

= -zSo sgnfi , (2.11) 

fi=IJ,=0 



where we have used the standard convention (dictated by the Vafa-Witten theorem) that 
the fermion condensate {0\ipip\0) and the corresponding quark mass are of the same sign. 
In this way the mass dependence of the fully quenched QCD3 partition function in the 
microscopic limit may be explicitly computed to be 

Zi'^) (/x, Ji)= e-' ^^« ('^-^^ ^^° ^ . (2.12) 



Later on we shall see that the simple form ( |2.12| ) agrees with what one obtains from the 



generic finite-volume field theory and also from random matrix theory arguments. 



3 Low Energy Effective Field Theory 

In the previous section we presented arguments in favour of a rather intricate pattern of 
fiavour supersymmetry breaking in partially quenched QCD3. This symmetry breaking 
pattern is crucial for the effective field theory computation of the microscopic spectral 
density. In the four dimensional supersymmetric formulation [0, the symmetry breaking 
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pattern is assumed to mimic as closely as possible the known bosonic one, and this turns 
out to be the correct answer. To help clarify the origin of the required finite volume 
partition functions, in this section we will derive the low-energy effective field theory for 
partially quenched QCD3 in the microscopic domain using the supersymmetry technique 
of random matrix theory |^. Although the main subject of this paper is to exam- 
ine properties of the Dirac operator spectrum using field theory methods, it will prove 
instructive to use random matrix theory as an intermediate step. This will provide, fol- 



lowing ^5], a rigorous derivation of the symmetry breaking arguments of the previous 
section. For instance, it will clarify how to choose the appropriate integration domain, for 
the definition of the partition function, as a Riemannian submanifold of the Goldstone 



supermanifold ( p..l6| ) 



For this, we start by writing down a random matrix model with the same global 
symmetries as the field theory (|1.12| ) in the Gaussian unitary ensemble. The partition 



function is 



^iv,,^„({m„};{/..,/i.}) = J DT e-^-^^ I] ^;' _ ]^JI ndet(T-.m,) . (3.1 



u{N) 



However, the matrix model (|3.1|) is difficult to deal with directly because of the asymmetry 
between the valence and sea sectors. It will prove convenient throughout this paper to 
have a completely supersymmetric form of the field theory, even in the physical sector. 
To this end we introduce a set of bosonic sea quarks which are the superpartners of the 
dynamical fermions. They have masses rhi, i = 1, . . . , Nf, which by supersymmetry are 
distributed according to sign as in ( p..l[ ) in the mass matrix 

TW = diag(mi,. . .,mn^,-mi,. . . ,-m„J . (3.2) 

The matrix model partition function of this extended supersymmetric field theory is 
defined byB 

^^ det(T - ifii) i^ det(T - imj) 



Z.,.„({m.,m.};{^.,/..}) = J DT e"^--^ n|^^^ R 



.=1 det(T - ifii) fj-^ det(T - irrij) ' 

(3.3) 

and the desired partition function (|3.1|) may then be computed from the infrared regime 
of the bosonic sea sector of ( |3.3| ), 

^A^/,Af«({"^i};{/"i)Aii}) = n -1™ {-ifhjfZj^j^i{mi,fni}\{\ii,lii}\. (3.4) 

We will denote by A^t = Nf + N^ the total number of flavours. 



^In the following the notation Z is strictly used for supersymmetric partition functions only. 
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3.1 Supersymmetric Representation 



We will now apply a colour-flavour transformation |26] to rewrite the integration in (3.3) 
as an integral over the total flavour space of the valence and sea quarks. Let us introduce 
an (A^tI^t) complex supervector by 



* 



>■ 



(3.5) 



where ipi, i = 1, . . . , A^t, are complex fermionic variables transforming in the vector rep- 
resentation of U{N) which can be thought of as representing the valence and sea quarks, 
while 0Q,, a = 1, . . . , A^t, are complex bosonic variables, also in the vector representation 
of U{N), which can be thought of as representing the ghost valence and sea quarks. The 
determinants appearing in (|3.1| ) may then be exponentiated in the form 



^ det(T — ifii) jt det(T — irrij] 
i=i det(T - ifLi) ^Ji det(T - irhj] 



-irr 



D^ D^ q-^'^\^®Mnt,\Nt,)+^n®M)-^ 



(3.6) 



where 



M = diag (M,iji, 



. /^Af„ 



M,n 



1) 



. /^Af„ 



(3.7) 



is the mass matrix (ordered appropriately according to sign as in (|1 . 1| ) ) , and the super- 
vector measure is deflned by 



N Nt 
a=l a=l 



■'a,a ^'ra,a 



n 



d 



{3.t 



The symbol A in (|3l6| ) denotes the graded inner product of supervectors ( p.Sj) deflned by 
while the adjoint \E' of the supervector ( |3.5|) is deflned as 



l„,,-l„,,l„+,-l. 



(3.10) 



^ = ^t f] ^ith Q = diag (Iat^ 

The deflnition ( p.lO[ ) ensures that the integration in ( |3.6| ) is uniformly convergent. 

The integrations in both (p.3|) and (|3.6|) are uniformly convergent. Interchanging them 



to perform the Gaussian integral over the Hermitian matrix T, we get 



^Nf,N^ [{mu rrii}; {yUj, /ij} j 



-^\nf+n^ 



D"^ D^i exp 



('C^tI'Vt') 



^A (l 



TV 



M]-^ 



2NT.I 



STr 



'Y.a^i®^a -Ea0l®^a' 



Ea^l 



Ea0i 



(3.11) 
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Here STr denotes the supertrace on (A^tI^t) supermatrices defined by 



Aff Abf 
Afh Ah 



STr " 7 = trA;^- trA, , (3.12) 



where Af^ denotes the bosonic fermion-fermion block, A\,\^ the bosonic boson-boson block, 
and y4f,j, Ajf, the Grassniann boson-fermion blocks of the supermatrix. These blocks are 
all A^T X ^T matrices. The sums in ( p.ll|) run from 1 to A^. 



To do the vector integration in ( p.ll ), we rewrite the four vector interaction term in 



Gaussian form by using the generalized Hubbard-Stratonovich transformation 



^STr(^^M..®*,) = 



NT.t 



DA e 2^STrA^+J<JA(^l]v®Aj* (3 13) 

where -DA is the invariant Haar-Berezin measure on the Lie superalgebra gl{Ni[^\N'j-) . The 
elements of this superalgebra may be parametrized as 

where A G u{Nt) contains ordinary mesons made of quarks and antiquarks, x S'lid x ^^^ 
independent complex-valued Grassmann matrices representing fermionic mesons consist- 
ing of a ghost quark and an ordinary anti-quark, and the boson-boson block iX, which 
represents a meson constructed from two ghost quarks, parametrizes the Lie algebra of the 
non-compact group GL{Nt^,C)/U{Nt). Because of the supertraces, this compact/non- 
compact structure is required for convergence of the integration in ( |3.13 ). The normalized 



Berezin integration measure in (|3.13| ) may be written in terms of the Cartesian coordinate 
parametrization ( p.l4|) as 



DA = (-vr)-^^/^ n ^A,- n dK, ® R H ^ ^ • (3-15) 

i,j=l a,/3=l k=l a=l ^Afccr l^Xcrk 

However, the resulting integration over the supervector ^ only converges when the 
boson-boson block A of the supermatrix A takes values in a non-compact u{nf +n^ ,nf + 
n~) subalgebra of gl{N'j^,C) defined by the condition A"'" = fiAfi.l3 Then the integra- 
tions are all uniformly convergent, and so by integrating first over the supervector \1/ in 
( p.ll| , p.l3D we arrive at 



/iVS AT 

DA e-^^T'-^'[SDet(A-z^)] (3.16) 

gl{NT\NT) 

where the superdeterminant is defined by 

\Afb Abb J det Abb det [Abb -AfbAjJAbf) 



*See |25|| for a more precise description of the entire integration domain. 
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and it satisfies STr In A = In SDetA. The expression ( p.l6| ), which is exact at the level 
of the random matrix model, is the desired representation of the partition function as an 
integral over the flavour superspace. Note that the partition function ( |3.16| ), in the super- 
symmetric limit and in the limit of massless sea quarks, is invariant under the transforma- 
tions \1' I— i> f/ ^ which leave the bilinear form \&A\& invariant. Such transformations satisfy 
the condition W = Q U^^ Q and form the unitary supergroup [/(A^tI^/ + '^J, nf + n~). 
This illustrates how the matrix model captures the global symmetries of the original field 
theory. In the next subsection we shall derive the symmetry breaking pattern from this 
point of view. 

3.2 Local Scaling Limit 

Let us now consider the partition function ZN^N^{{mi,mi};{fLi,fii}) in the thermody- 
namic limit. Precisely, we want to study the superintegral ( p.l6| ) in the local scaling limit 
N ^ oo whereby the quark masses are rescaled by the mean level spacing p(0) = N'Eo/ii 
and held fixed. In this limit, the partition function (|3.16|) may be expanded according to 



^Nf,N, [{mi, TTii}; {/ij, fXiYj 
= f DA exp STr 



gl{NT\Nrr) 



^^° A^ + AT In A - -^ A-^A^, + 0{- 



2 So ' VA^ 



(3.18) 



where A^s = NTiqA4, and we have implicitly restricted the integration in ( |3.18D to in- 



vertible A (This will be valid within a saddle-point approximation to follow). In the large 
A^ limit, the integral (|3.18|) is dominated by the stationary points of the function 



J^(A) = STr ^A"-lnA . (3.19) 

The saddle-point equation is 

S^A-A~^ = 0. (3.20) 

The saddle-point equation (|3.2CI|) is invariant under G'L(A^tI^t) rotations of the su- 
permatrix A. We shall therefore solve it first for diagonal A = Aq, and then determine the 
adjoint orbit of Aq under the supergroup GL(Nt\Nt) which will produce the full solution 
superspace of (|3.2CI|) . For this, we set all the Grassmann variables to zero and consider 
diagonal fermion-fermion and boson-boson blocks (Recall from the previous section that 
this is in fact all that is required to deduce the pattern of flavour symmetry breaking). 
The saddle-points are then 

Ao = -^f (3.21) 
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where T is an (A^tI^t) diagonal supermatrix with T = Uat^iat^. In general, the matrix 
( p.21| ) does not lie in the integration domain required for uniform convergence of the inte- 
gration over the boson-boson variables. We will require that, via Cauchy's theorem, the 
integration contour can be analytically continued into the saddle-point manifold without 
crossing any of the poles of the function SDet^(A — ijii). The signs of the eigenval- 
ues of Ao are then uniquely determined by the supermatrix Q, such that the restriction 
Aq = i7 Ao fi is satisfied. In this way we find that analyticity and the forced choice of inte- 
gration domain for A fix the boson-boson part of T to be diag(l„^,, — In^, 1„+, ~^n')- '^^^ 
fermion-fermion block of F, for which there are no convergence nor analyticity constraints, 
depends on which configuration will dominate the superintegral in the limit N —^ oo. We 
will now determine this block using supersymmetry. 

For this, we expand the function ( p.l9| ) about the critical point (|3.21|) and evaluate 
the resulting Gaussian fluctuation integral in ( p. 18 ). We have 



^(Ao + Q)= ^(Ao) + ^ STr (g2 + f Q f g) + O [Q'') (3.22) 

where Q G (^/(A^tI^t)- Since the supermatrix F defines a projection operator on the linear 
space f7/(A^T|^T), we can make an orthogonal decomposition Q = Qc + Qo (with respect 
to the inner product {X,Y) = STrXY) corresponding to the ±1 eigenspaces of F, 

fQef = Qe , TQ,T = -Q,. (3.23) 

Then the quadratic fluctuations ( |3.22| ) depend only on the even degrees of freedom Qe, 

^(Ao + Q)= ^(Ao) + S2 STr Ql + O [Q^) . (3.24) 

The invariant Berezin measure factorizes as DA = DQc DQo. Upon substitution of (|3.24|) 
into ( p.l8| ), the integration over the Gaussian fluctuations Qg around the saddle-point 



produces one factor of A^~^ (resp. A^+^) for each commuting (resp. anticommuting) 
direction of steepest descent. The limit N -^ oo oi ( p.l^ ) is therefore dominated by those 
extremal hypersurfaces which have maximal transverse super-dimension d-^ — dj^. This 
dimension is zero when the fermion-fermion and boson-boson blocks of F are identical, 
and then the Gaussian fluctuation integral over the modes Qe produces unity in the limit 
N ^ oo. Therefore, the dominant saddle-point configuration is given by the supermatrix 

f = diag (l„^, -1„^, ]!„+, -]1„- |]1„^, -Ir^f, 1„+, -K- ) • (3-25) 

Having integrated out the degrees of freedom Qe transverse to the saddle-point super- 
manifold, let us now focus on the degrees of freedom Qo tangent to it. The stabilizer sub- 
group of the matrix ( |3.25| ), i.e. the group of matrices U G GL{Nic\Nt) with U T U~^ = F, 
is isomorphic to GLijif + n'^\nf + n^) x GLijif + n~\nf + n~), and so the saddle-point 
supermanifold, defined by the adjoint G'L(A^tI^t) orbits of ( |3.21| ), is the coset superspace 

G[nf- n„ , n J - c'L(ny + n+|n; + n+) x GL{nf + n-\nj + n") ' ^^'^^^ 
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The measure DQo is the local invariant Berezin measure of this space at the point Aq. We 
now substitute the adjoint orbits A = f/ Aq U~^ into ( p.l8| ) and use the fact pO|] that the 
integration measure DK coincides with the invariant measure for integration over the coset 
space ( [3.261 ). By using J^(Ao) = 0, we then arrive at the expression for the matrix model 
partition function (|3.3|) in the microscopic limit which has a hyperbolic supersymmetry, 



Z^^^^X{m,,m,}■{^^,,J^,}) = j DU e^^^o STr(>(^f c/-) (327) 

where here DU denotes a normalized invariant Haar-Berezin measure on the Lie super- 
group GL{Nt\Nt^). 

Finally, to obtain the partition function for the partially quenched QCD3 partition 
function in the microscopic domain, we need to take the large mass limit ( |3.4| ). This limit 
is a little subtle and is described in appendix A, where it is shown that (|3.27|) reduces to 



^J^k(i"^^}'i/^i'/^*} 



DU Q-'^^^o^'Tr{M'Ut'U-'') 



(3.28) 



where 



jCi' = diag (mi 



g{nf;nZ,n^ ) 



''TT'Nf,fJ'l, ■ ■ ■ ,/^Af„ 



/^l, 



<^^N^ 



(3.29) 



(3.30) 



is the mass matrix of the quantum field theory ( 1.12|) , and 

f = diag (l„,, -!„,, 1„+, -]1„- |l„+, -1„ 

Heuristically, the reduction of ( |3.27|) to (|3.28| ) in the regime where the modes of the 
ghost sea quarks are irrelevant can be understood by appealing once again to the Vafa- 
Witten theorem. The Hubbard-Stratonovich variable A which was introduced in ( |3.13|) 
may be thought of as the corresponding supermatrix of quark condensates, and there- 
fore, under the assumption of spontaneous flavour symmetry breaking, it should have 
eigenvalues of the same signs as those of the corresponding mass matrix ( p.7| ). Carrying 
out the colour-flavour transformation directly for the partition function ( p.lD itself (see 
appendix A), a repeat of the saddle-point analysis of this subsection using these symme- 
try breaking arguments to determine the dominant configurations would lead exactly to 
the effective field theory (|3.28|) . Indeed, the integration domain ( p..l6|) in (|3.28|) is the 
Goldstone manifold for the flavour supersymmetry breaking pattern (|1.15| ) , and the finite 
volume partition function ( p.28| ) is the appropriate generalization of the Itzykson-Zuber 
type integral ( |1.6| ). Note that the integration domain for the low-momentum Goldstone 
modes fits into the Zirnbauer classification of the local scaling limits of random matrix 
theories |^. This Riemannian symmetric superspace is supported by the compact sym- 
metric space U{Nf + N„)/U{nf + n^) x U{nf +n~) in the fermion-fermion sector, and by 
the non-compact symmetric space U{ri^ ,n':^)lU{ri^) x U{n~) in the boson-boson sector. 
These integration manifolds are defined by the intersection of the adjoint orbits of Aq 
above with the forced integration domain for A. 
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4 Quenched Approximation 



As a warm-up to the general case, in this section we will study the finite volume partition 
function in the fully quenched limit Nf = 0. Since our ultimate goal is to obtain explicit 
expressions for the microscopic spectral density ps, we will concentrate on the case of 
only a single species of valence quarks, Ny = 1. This particular case can be worked out 
in complete detail and it will serve to illustrate some of the general formalism of the 
previous section. In particular, it will shed light on the role of supersymmetry in the 
various manipulations. For example, a crucial issue within the present formalism is the 
supersymmetric limit of the effective field theory ( |3.28|) . According to the general setup, 
in that case the partition function of partially quenched QCD3 should reduce to that of 
ordinary QCD3, 



^J^,k(i"^^};i/^^'^^}) 



{Mi=Mi} 



7LS 



;(-M) 



(4.1) 



This reduction of the supersymmetric integral (|3.28|) is highly non-trivial and is related to 
the notorious boundary ambiguities, or Efetov-Wegner terms [^ , which plague integrals 
over super-manifolds. They are related to non-integrable singularities of the Berezin mea- 
sure at the boundary of the integration domain in certain coordinate parametrizations 
whereby the volume form of a non-compact supermanifold becomes a form- valued differ- 
ential operator. Some formal mathematical descriptions of these boundary terms can be 
found in [25, 28, 29 1. In appendix A we demonstrate the reduction (^4.1|) formally starting 



from the finite A^ representation (|3.16|) . In this section we shall see how it comes about 
through explicit calculations. 

We first consider the saddle-point analysis of the previous section in this special in- 
stance. There we saw that the dominant configuration in the large A^ limit is determined 
by the supermatrix T = sgn(yu) li|i, for which Q = Q^. The main consequence of this 
result is that the saddle-point supermanifold is a single-point, since U KqU~^ = Aq for 
all U G GL(1|1), and the partition function localizes onto its value at the unique critical 
point Aq = ^TT-^ li|i. By substituting this solution into (p.l8|), we arrive at the final result 



for the fully quenched QCD3 finite volume partition function in the static limit. 



Zoj\f^,f^) = e 



-isgn{/i) STrMs 



(4.2) 



which coincides with the result ( p. 12 ) obtained directly from the quantum field theory. 



(oo). 



Note that in the degenerate case /i = /i, we have Zqi\p,p) = 1, as expected. These 
results simply refiect the fact that, in the microscopic regime of fully quenched QCD3, 
there is no spontaneous breakdown of the G'L(lll) fiavour supersymmetry. 

Let us now go back and consider the partition function ( p.l6| ) in the quenched limit. 



Zo,iip,p) 



DA 



-^STVA^ SDet^ 



glim 



A 




(4.3) 
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Note that for fi = fi, the integrand of ( |4.3| ) is a supersymmetric invariant function, so 
that the Efetov-Wegner theorem imphes Zo,i(/i, /x) = 1 (see appendix A, eq. ( |A.10| )). On 
substituting in the parametrization ( |3.14| ) and integrating over the Grassmann variables, 
eq. (^73|) becomes 



^0,l(/i, /i) 



N 



jN 



n 



d\ 



dX e 



A^+A^ 



(A — i/i) 



N-l 



y\- fi + i sgn(/i) e 



N+l 



-rill 



(A — in) 



N 



[X- n + i sgn(/i) e 



N 



(4.4) 



where the parameter e — > 0"^ regulates the poles of the integrand at A = /2. Note that, 
according to the general analysis of the previous section, the particular choice of analytic 
continuation of the integration domain into either the upper or lower complex half-plane 
depends on the sign of the valence masses. This sign dependence is required for conver- 
gence of the Hubbard- Stratonovich transformation of section 3.1. 

The integrations over A and A in ( |4.4| ) decouple. The A integrals can be evaluated in 
terms of Hermite poljTiomials Hn[x) BO] by using the integral representation 



(2z) 



H^{x) = ^L / dt {t-ixY e 



n > 



(4.5) 



For the A integrals, it is convenient to rewrite them in terms of shifted Gaussian moment 
integrals by using the Fourier transformation 



1 



\ — fi + i sgn(|U) e 



(isgn/i)" 

{n-iy. 



dk k^~^ ^isgn{ii) k[\-fi+i sgn 



M^) 



(4.6) 



which is valid for n > 0. By substituting ( [4.6|) into ( [4. 41) , and performing the Gaussian 
integrals over A, the remaining integrations over k can be expressed in terms of the 



generalized Hermite functions ||31 



TT 



dtf" e 



-r-2Jxi 



n> 



(4.7) 



The Hermite functions are non-polynomial and they are related to the error function. 
Their imaginary parts coincide with the Hermite polynomials ( [4.5|) , while their real parts 
represent the second set of linearly independent solutions of the Hermite differential equa- 
tion which can be expressed in terms of confluent hypergeometric functions as 



ReTCJx) 



(-l)'+Si^iU-fc;i 



3.2 

1 2 ' 



' 2 ' 



X 



n 



n 



2k 

2k -1 . 



{U 
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Combining these results, we arrive at the exact expression for the finite volume partition 
function of fully quenched QCD3 with one species of valence quarks, 



Zo,iil^,fx) 



TT 



2^{N-1] 



-iVS2/i2/2 



H 



N 



- H 



N-l 



Afs: 



^fx\n 



N 



N^ 






/^l 



NY^ 



(4.9) 



Let us now demonstrate that ( [4 .91 ) leads to the properties of the finite volume partition 
function deduced by formal arguments above. First of all, in the degenerate limit fi = fi, 
we may use the generalized Christoffel-Darboux formula to deduce i3l[] 



Hn{x)Hn-l{x) - nn{x)Hn-l{x) 



[n 



1)! 



TT 



(4.10) 



which when applied to ( |4.9D leads to Zo,i(/i, /i) = 1, as anticipated. This demonstrates 
once again that the valence mass independence of the degenerate partition function is a 
highly non-trivial result of the exotic properties of superintegrals. Secondly, let us check 
the microscopic limit of the partition function ([4.9| ). For this, we need the asymptotic 



form of the Hermite functions 31 



TiJx) 



— > en 



TT n—\J1nix 



and the standard asymptotic forms of the Hermite polynomials |3^ 

(-1)'=2'=(2A;-1)!! e^'' cos (ViA^TT 



EJx) 



4)^'2'=-5(2A;-2)!!V2A;-l e^' sin (V4A; - 1 



X 



(4.11) 



n = 2k 

n = 2k~l . 
(4.12) 



We substitute ( 4.11| ) and ( 4.12| ) into (O) and take the N —>■ 00 limit with the rescaled 
masses NT^q fi, NT,q ft fixed. By using the Stirling approximations 



n\ 



72 



nn rf e " 



n! 



27r n"/2 e -"/2 ^ 



(4.13) 



we then find that (|4.9| ) in the local scaling limit reduces exactly to the anticipated result 
( p.l2| ) for the microscopic partition function in this case. 



5 Microscopic Spectral Density 



We finally come to the evaluation of the microscopic spectral density ps{u;u!i, . . . ,ujnj,) 
from the finite volume, supersymmetric field theories (|3.28|) . We will compare the expres- 



sions obtained from this analysis with those computed in [|1^, |14[ using random matrix 
theory techniques. In the next section we shall generalize this analysis to compute all 
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microscopic /c-point spectral correlation functions. This will thereby demonstrate that 
the supersymmetric formulation of partially quenched effective field theories provides an 
analytical framework in which one can establish the equivalence between the microscopic 
Dirac operator spectrum of QCD3 and the microscopic spectral correlators of random ma- 
trix theory. We shall begin with the quenched approximation to QCD3, and then move 
on to the general case of Nf > flavours of dynamical fermions. 

5.1 Quenched Limit 

In the quenched approximation Nf = the partition function in the microscopic domain 
is given by ( |2.12D . The valence quark mass dependence of the fermion condensate may be 
computed by using (|1.13|) to get 



S,(n) = So--. (5.1) 



The function ( |5.1| ) has a jump discontinuity of 2i Sq across the real axis. By using ( |1.14| ) 
we thereby find that the microscopic spectral density is given by 

Ps{u) = - , (5.2) 



which is the expected result in this case |TD|, |n[. The spectral distribution function 
is fiat and it coincides with the usual macroscopic density p(A) evaluated at the spectral 
origin. Of course this result is not unexpected, given the absence of dynamical fermions. 
The eigenvalue density near the zero mass regime would normally contain an oscillatory 
fine structure with period set by the mean level spacing vr/iVEo. In the present case, this 
fine structure is absent, and the eigenvalues of the fully quenched QCD3 Dirac operator 
are on average uniformly distributed over the real line. 

5.2 QCD3 with Nf Flavours 

To treat the general case, we need an appropriate parametrization of the Goldstone man- 
ifold Q{nf; l'*^'^^,0) in ( |1.16|) . The ordinary integration manifold supporting this coset is 



the symmetric space U{2nf + 1)/U{nf) x U{nf + 1) in the fermion-fermion sector, while 
it is simply a point in the boson-boson sector. Unitarity then requires that the (2n/ + 1)- 
dimensional anticommuting vectors x ci-nd x which comprise the Grassmann components 
of the corresponding Goldstone superfields U obey the constraint xU~^ x = 0, where 
U are the commuting, unitary fermion-fermion degrees of freedom of U. Because of the 
asymmetrical decomposition of this vacuum manifold, it is difficult to evaluate the integral 
( ^.2^ ) in a straightforward way. 

To circumvent this asymmetry, we proceed in a way that is reminescent of the observa- 



tion 1^ that the microscopic spectral density is related to the finite volume effective field 
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theory for QCD3 involving two additional species of quarks of equal imaginary mass. By 
using ( |1 .131 , 11 -I'D it is straightforward to show that the spectral density in the mesoscopic 
domain of the underlying finite volume gauge theory may be computed directly from the 
partially quenched partition function with two flavours of valence quarks as 

Ps{u;uJi,. . .,UNf) 

1 1 .. e d d ^ f M, 



e^0+ TT Olii Olio VA/En / 



iVSo Z;v,,o (^) ^-0+ ^ dfi, dfi2 ''"' ViVEo' 



(5.3) 

In ( |5.3|) the masses /ij, fti are all positive initially and then analytically continued into the 
right complex half-plane. To prove the identity (5^), we note that by using ( 1.12| ) the 
right-hand side may be computed to be 

_ lim i V i (- 1 1 ) (5.4) 

^^0+ TT ^ Am - A„ - 2ze \Xn- X + te Am- X-teJ 

where A„ are the Euclidean Dirac operator eigenvalues which are assumed to be non- 
degenerate. Because of the overall factor of e in ( |5.4| ), the m ^ n terms each vanish in 
the limit e — * 0+. This same factor is cancelled by each of the m = n terms in ( ^.4| ), 
which when summed reproduce the expression ( |1 . 14|) . Again it is convenient to evaluate 
the partially quenched partition function in ( ^.3|) by using a completely supersymmetric 
expression for it. As in section 3, this is achieved by introducing very heavy superpartners 
to the dynamical fermions with mass matrix ( |3.2|) . Given that the large mass and N ^ 00 
limits commute (see appendix A), to treat the thermodynamic limit it is convenient to 
write the large mass expansion in terms of a ratio of finite volume partition functions as 



M' 



zi^l{M,M) 



The partition functions appearing on the right-hand side of (|5.5| ) can now be readily 
analysed. 

Coset Parametrization 



The integral ^jv^ 2 is given by (|3.27| ) defined over the symmetric superspace G{nf\ 1, 1) 



in (|3.26| ). To parametrize the supermatrices U of this coset space, it is convenient to 



change basis on the underlying supervector space C^^Ia^t ^ (^n^+iin^+i ^(^n^+iin^+i ^^ ^^ 



orthogonal decomposition into the ±1 eigenspaces of the projection operator (|3.25| ). In 
this basis, the parity matrix T takes the form 

f = l„^.+i|„^.+i ® as (5.6) 
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while the mass matrix 



with 



becomes 

M =diag(Mi,-M2) 



(5.7) 



' Snj+1 Si ) 






mi,. . . ,mnf.,fii 



{5.t 



M, = diag(cf\ '^'^ 

= diag (mi,..., m„^,/ii 

for i = 1,2. By exponentiating the coset generators it is then straightforward to show 
that the matrices U G G{nf] 1, 1) can be parametrized as 

/Vl + TT T \ 

~ l^ T Vl + TTJ ' 

where the supermatrices T, T G GL{nf + l\nf + 1) are related by 



(5.9) 



T = diag(lny,+i 



ln,+l) Tt . 



(5.10) 



In these Cartesian coordinates, the invariant measure for integration over the coset space 
is given by 



ny+1 nf-\-l Ytj-T-L nj-]-L o 

Du= II rfT,, dT*^ n ^T^-/^ ^T^:/. ® n n 

i,j=l a, 13=1 



d d 



d 



\ ^A dTk. 5TL dT^k 5T:, 



-. (5.11) 



However, despite the simplicity of the integration measure (|5.11|) , Cartesian coordi 
nates are not convenient for the evaluation of the integrals ( p.27|) . Following |29, 3^ 



we introduce the Efetov polar coordinate parametrization of the coset space ^^. These 
coordinates are inherited from the decomposition of [/(A^tI^t) matrices into eigenvalue 
and angular degrees of freedom. When projected onto the coset, these variables form an 
orthogonal decomposition into parity even and odd sectors. Namely, we may parametrize 
the elements of the coset space as 



U = VEV 



-1 



(5.12) 



where the anticommuting coordinates reside in the angular, "eigenvector" matrix V which 
commutes with the coset projection operator. 



VT = TV 



(5.13) 



while the commuting coordinates live in the "eigenvalue" matrix S which anticommutes 
with the parity matrix. 



sr 



"Z'-l 



(5.14) 



From ( p.l3| ) it follows that the angular matrices V G U{nf + l\nf + 1) x U{nf + l\nf + 
1)/U{l)^f~^^^"'f~^^ admit the matrix presentation 



y = diag(V^+,1/_ 



(5.15) 



23 



in the parity ordered basis introduced above. The matrices S satisfying (|5.14|) may be 
parametrized as 



in+i 



ITZ-l 



7^-l 



7^+l 



where 



7^ 



diag( 



n 



} ■ ■ ■ 1 'rif+l 



ri 



} ■ ■ ■ 1 lUf + l 



(5.16) 



(5.17) 



The compact fermion-fermion radial coordinates r^ each hve in the finite interval [—1, 1], 
while the non-compact boson-boson radial coordinates f^ live in the semi-infinite interval 
[l,oo). The collection of matrices ( |5.16| , |^TT7| ) form an [uf + l\nf + 1) maximal abelian 
subgroup for the Cartan decomposition of GL{Nt\Nt) with respect to the stability sub- 
group GL{nf + l\nf -|- 1) x GL{nf + l\nf + 1). We recall from section 3.2 that the latter 
degrees of freedom defined the directions of steepest descent on the saddle point manifold, 
while the former ones determined its structure. 

There are two main advantages of this polar coordinate parametrization. First of all, 
while the original integration in (|3.27|) cannot be trivially extended from the Goldstone 
supermanifold to the full fiavour supergroup, the angular integrations over V± can be 
extended to the whole unitary supergroup U{nf + l|n/ -|- 1). Then, the appropriate 
supersymmetric generalization of the Itzykson-Zuber formula may be applied. Secondly, 
the choice of coordinates ( p.l2| ) will completely decouple the bosonic and fermionic sectors 
of the supergroup integral from one another in such a way that the large mass limit ( |5.5| ) 
may be easily taken. The price to pay for the introduction of these coordinates is that 
the radial integration domain has a boundary, and so we can anticipate the appearence of 
Efetov-Wegner terms. The calculation of the Berezinian of the coordinate transformation 



( ^.12| ) can be done in the usual way [^ and the measure assumes the familiar form 






{5.U 



i=l 



«=1 



where DV± are invariant Haar-Berezin measures on U{nf + l\nf + 1) and the A's are 
{uf + 1) X [nf + 1) Vandermonde determinants which are defined in appendix B. Here and 
in the following we will, for simplicity, not keep track of numerical integration factors and 
simply denote them collectively by Co- The appropriate normalization will be restored 



by hand in our final result later on. From ( |5.18|) we see that the Berezin measure in 
polar coordinates contains non-integrable singularities at the points r^ = f^ = 1 for 
any pair (i,a), unlike the analytic Cartesian coordinate measure ( p. 11 ). These fictitious 
singularities are caused by the mixing of nilpotent terms into the commuting degrees of 
freedom in ( |5.12D which lead to total derivatives that give rise to additional contributions 
to the pertinent integral at the boundaries of the radial integration domain. Although 
these boundary terms can be calculated in principle by using the techniques described in 
29| , they are rather cumbersome in form and not very informative. In what follows we 
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will show that, as in ||33|, they do not contribute in the limits of interest. Heuristically, 
these terms arise from the introduction of the fictitious superpartners to the physical 
quarks and the valence fermions. They therefore should not contribute to any physical 
observable, such as the Dirac operator spectrum. 

Finite Volume Partition Functions 

We will now evaluate the partition function Z^ 2 t'Y simply ignoring the Efetov-Wegner 
terms. We substitute (|5.12|) , (|5.15| ), ( ^.16| ), and ( ^.18| ) into (|3.27| ). By using the commu- 
tation relations ( p.l3| ) and ( |5.14|) we find that the supertrace factorizes into parity sectors 
as 



STr (MV^V-^YV^-^V~^') = STr (^MiV+nV'^) + STr (M2F_7^y^l 
and so the partition function can be expressed in a factorized form as 



(5.19) 



^i^,2 (A<,A^;yUl,-/i2,yUl,-Ai2 



«/+! I 



rif + l 00 



X f DV+ e -' ^^« ^'^' (*^^ ^+ ^ ^+ ') 

U{nf+l\nf+l) 

X / DV- Q-i^^oSTT{M2V-nvz^) 

U(nf+l\nf+l) 



(5.20) 



The unitary integrals in (|5.2CI| ) can each be evaluated by using the supersymmetric gen- 
eralization of the Itzykson-Zuber formula ||2^ (see appendix B, eq. ( |B.13| )) to get 



-2j^,2 (M,M; Hi, -H2, /^i, -/X2) = C'o n 



A 



■^(fc)^^(fc) 



fc=l,2 



A 



((fe) A C^'') 



rif+l 1 rif+1 00 



i=l 



1 "=1 1 



fc'=l,2 



«J 



-JAfEoCf'V, 



det 

a,/3 



,iNEoC2''^rp 



(5.21) 



The radial integrals over r^ and r^ in ( 5.21 ) decouple. Because of the permutation symme- 
try of the integration measure, the two fermion-fermion determinants may be combined 
into a single one detij e~* °^^* "'"''j ''^^ , times the order (n/ + 1)! of the permutation 
group Snf+i which we absorb as always into the normalization constant Cq. The same is 
true of the boson-boson determinants. The radial integrals may now be straightforwardly 
done, and by using the definitions of Q and Co in (|5.8| ) we arrive at 



Z 



Nf,2 



{M,M; 



/^l) ~/^2)/^l; ~f^2 
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Co (/ii-/xi)(/x2-/i2) n 



^ {rrii - /ii)(mi - 7x2) -rC (m„ - /ii)(ma - /i2) 



i=l 



(rrii - yUi)(mi - /ia) iJi (m« - /ii)(m„ - /i2) 



X 



A[m]2A[m]2 



det 



sin A^Sq (rrij + mj) sin A^Sq (mj + fii) 

NT,o {mi + /ii) 

sin A^So (/ii + /i2) 

A^So (/ii + /i2) 



X det 



NTjq {rrii + rrij ) 

sinA^So (/i2 + "^j) 

A^So (/i2 + "^i) 



A^So(mQ+m/3) A^So(ma + /ii^ 



, JiVEo (Ai2+»n.^) 



, i7VSo{Ail+A2) 



A^So (/i2 + "^/3) A^So (/il + /i2) 



(5.22) 



In (|5.22|) the A's denote Uf x Uf Vandermonde determinants in the fermion masses rrii and 
ffia, i,a = 1, . . . , Uf, while the ordinary determinants have dimension {nf + 1) x {nf + 1). 

Let us now evaluate the partition function 2j^q which appears in the denominator 
of the expression ( ^.5| ). Applying the exact same steps which led to ( p.22| ) for the coset 
integral (|3.27|) over G{nf] 0, 0), we arrive at 



Z 



(oc) 
Nffl 



M,M 



Co 



A[m, m 



^12 



det 
A[m]2A[m]2 i<i,j<nf 



sin A^Eq [mi + rrij 



NHq {rrii + fTij) 



X det 

l<a,/3<n^ 



A^So (m„ + rhp) 



(5.23) 



The terms in ( |5.23|) which depend only on the physical fermion masses rrii are readily rec- 



(00) 



ognized as the finite volume partition function Z^'q{M.) for ordinary QCD3. Indeed, one 
can parametrize the ordinary symmetric space ( |1 . 7| ) in the same manner described above 
and evaluate the finite volume partition function ( |1.6D analogously as an integral over the 
coset space, rather than the full unitary flavour symmetry group, by using the ordinary 
Itzykson-Zuber formula for U{nf) |3^. In this way one may arrive at the representation 
(up to an irrelevant normalization factor) 



Zm,{uji, 



'Nf 



,^Nf) = ^ r no det 

■^^ A[CJ]2 l<i,j<nf 



sinh(ci;j + uJj) 

Ui + UJj 



(5.24) 



where here we have analytically continued the unfolded masses (|1.8|) to imaginary values 
to facilitate comparison with previous results. In [1^ an expression for the finite volume 
QCD3 partition function was derived by applying the ordinary Itzykson-Zuber formula for 
U{2nf) in a suitable regulated limit that removes the n/-fold degeneracy of the eigenvalues 
of the matrix T^. The expression ( p.24[ ), which utilizes the same integration formula but 
does not require dealing with any degeneracies, is much simpler and compact as it involves 
elementary rif x n/ determinants, rather than the 2n/ x 2n/ determinants that appear in 
|10| . While we have no direct proof at present that these two expressions are equivalent, 
we have checked that they agree in a number of cases. The results ( |5.23|) and ( |5.24|) clearly 
show that the Efetov-Wegner boundary terms vanish in the decoupling limits m^ — i> 00, as 
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expected. For example, one of these boundary terms comes from evaluating the integrand 
of ( |3.27|) at the origin U = l.^^\j^^ of the coset superspace ||29[. This adds the term 
g-jAfSo STr(Mi+M2) |.q j.^^^ above expressions, which produces a vanishing result in the 
large mass limit. 

We can now finally write down the desired expression for the partially quenched finite 
volume partition function. By substituting ( ^.22| )-( p.24| ) into ( |5.5| ) and taking the limit 



mn 



oo, we arrive at 



,(oo) 



^7V;,2(-^;/^b-/^2,/il,-Ai2) 

(/il-/il)(/i2-/i2) e^^^°(^i+^2) 



Co 



X det 



A [ml 



n 



{rrii - fii){mi - ^2) 



A^So (/ii + Ai2) r=i (mi - iii){mi - /^a) 



sin A^Sq {rrii + rrij) sin NT,q {rUi + /Xi) 



A^So {mi + nij) 

siniVSo (yU2 + ^3] 

NT^o (/i2 + rrij) 



A^So(mi + /ii) 

sinA^Eo (/ii + /i2) 

A^Eo (/ii + /i2) 



(5.25) 



The expression ( p.25| ) of course represents only the bulk, regular contribution to the su- 
persymmetric integral (p.27|) . The anomalous boundary terms are obtained by setting 
"f^i = 'f'a = ^ for one or several pairs of radial coordinates {ri,fa)- In the original integral 
( p.27| ), this corresponds to setting some of the supersymmetric blocks of the unitary ma- 
trix U equal to the identity matrix of the appropriate dimensionality. The integrand of 
( p. 27] ) then becomes a supersymmetric invariant function of the given block, and the inte- 
gral becomes correspondingly dimensionally reduced (see appendix A). In the integration 
measure ( p.l8| ) the Vandermonde determinants are reduced accordingly by omitting the 
given singular factors. The reduced coset integral can thereby be evaluated as above. As 
we have shown, only the boundary terms which are associated with the valence fermions 
contribute. These terms can be computed using the results of [^|. One of them comes 
from setting the (2|2) valence block of U equal to the identity matrix. Following the 
derivations above, this produces the boundary term 



^iVf,2(-^;/^b-/^2,/Xi,-/i2) 



A [ml 



-iNSoifii 



-Ml) gJA^So (^42-/^2) 



X det 



sin A^Sq (mj + rrij) 
A^Sq {rrii + rrij) 



which is responsible for the normalization (|4.1| 



(5.26) 



Spectral Distribution Function 

Finally, we can now easily compute the spectral density of the QCD3 Dirac operator in 
the mesoscopic region using ( ^.3|) . We note first of all that only the regular part (|5.25|) 
contributes to the spectral density. The factor of (/ii +^2)'^ that appears there is needed 
to cancel the factor of e in ( |5.3|) in the limit fii + fi2 —>■ 0. The boundary terms contain 
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fewer non-compact integrations, so they either vanish in the hmit ffia — > cxo, or they do 
not contain this factor and so vanish in the hmit e ^ O"*". With this in mind, we can 
now differentiate the expression ( p.25| ) with respect to the valence quark masses /Xj and 
take the hmit dictated in ( |5.iJ| ). By using ( |5.24D , we then arrive at a relatively simple 
expression for the microscopic spectral density, 

sinh(a;j + uJj) sinh(a;i + iu) 



det 



Ps[U\Ui, 



,^Nf 



UJi + UJj 



sinh(c(jj 



lU] 



UJi + iu 
1 



UJj 



lU 



TX 



det 



sinh(ti;j + ujj) 

UJi + UJj 



(5.27) 



where we have again analytically continued both the unfolded masses uji and the unfolded 
Dirac operator eigenvalues u to imaginary values, as these are the standard conventions 
that are used in gauge theory computations. The determinant in the numerator of ( |5.27| ) 
is of dimension (n/ + 1) x {jif + 1), and the overall normalization constant Co = l/vr 
has been fixed by the usual matching condition between the microscopic density and the 
macroscopic density at the spectral origin, 

P(0) _ 1 



lim ps{u]Ui,...,ujNf) 



(5.28) 



5.3 Examples 



The universal double-microscopic spectral density was computed in [T^ using random 
matrix theory techniques and found to be given by a rather involved determinant for- 
mula involving Bessel functions of half-integer order. Here we have found that the su- 
persymmetric method based entirely on the field theory formulation leads to an elegant 
expression (|5.27|) for the same quantity which is much more compact and convenient to 
use. Again, we have no direct proof of the equivalence of these two representations of the 
spectral density, but we note that they both involve determinants of the same dimension 
(jif -|- 1) X {rif + 1). We have checked that they agree in a number of special cases. For 
example, consider the case of two physical quarks, nj = 1, of equal and opposite mass m. 
By using the trigonometric identity 



2 sinh(x -|- y) sinh(x — y) = cosh 2x — cosh 2y 
the resulting 2x2 determinant in ( ^.27|) can be worked out to give the result 



Ps{u;uj,-uj) = - 1 + 



uj cos 2u — cosh 2uj 



n 



u^ 



uj'^ sinh 2uj 



(5.29) 



(5.30) 



which agrees with the known density of states from random matrix theory [0. Similar 
computations can be done for higher numbers of massive fermion flavours, and in each 



case we have found precise agreement with the results of |10 
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An important special case that can be worked out straightforwardly is that of an ar- 
bitrary number Nf = 2n/ of massless quarks. In this case, the ratio of determinants in 
( p.27| ) produces an indeterminant form and must be defined by an appropriate regulariza- 
tion in the limit Ui —>■ 0. In that limit, the argument of the determinant in ( [5.24] ) admits 
the Taylor series expansion 



f{uJi + UJj) 



= 2^ Aki Bii Bkj + 0{uj,' 



UJi + Uj 



kA=l 



(5.31) 



where we have defined the Uf x Uf matrices 



•^j J 



(i+i-2) 



(0) 



B. 



uj: 



i-i 



u 



(i-l) 



(5.32) 



In writing ( p. 31 ) we have used the fact that both the numerator and denominator of ( 5.24| ) 
vanish as a;^^ in the limit Ui — > 0. Since deti3 = A[u;]/nrii(* ~ 1)'; ^^is regulates the 
partition function (|5.24|) and yields the finite result ^J^^(0, . . . , 0) = deiA/{Y]^li{i-l)\Y. 

By substituting ( ^.31|) into ( ^.27|) , the microscopic spectral density can be written as 
the [uf + 1) X {rif + 1) determinant 



p,(u;0,...,0) = - det ( ^^"^^ 
where we have defined the n /--dimensional vector 



B-^Aa 



r 



f 



(i-i)i 



zu) 



(5.33) 



(5.34) 



The determinant ( 5.33| ) is readily evaluated by performing minor expansions along the 
last row and column with the finite result 



Ps(u; 0,...,0) = -(l-a^ Aa) . 

7]- V / 



(5.35) 



By substituting ( |0^ ) and ( ^341) into (|05| ), and by using /(2"-i)(0) = and /(2")(0) 
l/{2n + 1), the spectral density in the massless limit can thereby be written as 



Ps{u; 0,...,0) = - 



n 



Uf-l 



1-E 



1=0 






d 



k-l 



siiau 



d' 



'2l~k+l 



smu 



21 + 



fc=i 



du^ 



u 



du^l-k+l 



u 



(5.36) 



The function ( |5.3(j| ) can be expressed in terms of regular Bessel functions Jy{x) of half- 
integer order i/ [30|| by using the derivative formula 



4+iW 



2r^""M-h3 



2 / d 

-K \ dx^ 



smx 



X 



(5.37) 



which is valid for positive integral n, and the three-term recursion relations 



2u 



X 



2Jl{x) = Jy^i{x) - Ju+i{x) , 
Jy{x) = J^_i(a;) + J^+i(x) . 



(5.38) 
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After some algebra, the expression (|5.36|) can be simplified to the compact form 

u 
4 



U V 
Ps{u; 0, . . . , 0) = - [J„^_i (m)2 - Jnf + liu) Juf-liu) + Jnf + liuf - 4^_l («) 4^ + |(m) 

(5.39) 



which coincides with the massless spectral density obtained originally from random matrix 
theory [^, |36[. Here we have derived it directly from the effective finite volume partition 
function of QCD3 in the microscopic scaling regime. 

In [|T^ it was shown, by matching exact results from random matrix theory with the 
low energy effective field theory (|1.6| ), that it is possible to express the spectral density 
of the QCD3 Dirac operator in terms of a ratio of two finite volume partition functions, 
one of which involves two additional fermion species of equal imaginary mass, as 

Ps{u-u;i,...,u;nj) = —[[[u +u;,) ' ls/, , — ^ • (5-40) 

As mentioned at the beginning of section 5.2, this is not unexpected as the spectral 
density is given directly from ( ^.3|) . The partition function in the numerator of ( |5.40|) 
is understood as an analytical continuation in the additional fictitious fermion masses 
in which both mass parities are substituted by the value iu. The expression ( |5.27| ) can 



thereby be thought of as an explicit realization of this feature, with the appropriate 
analytic continuation of the coset representation ( |5.24| ) given by the numerator function in 



( p.4(j| ). This is the biggest advantage of the polar coordinate parametrization of the coset. 
The supersymmetric partners to both the sea and valence quarks are just spectators in this 
formalism, and it is the valence fermions themselves which give the explicit representation 
of the microscopic Dirac operator spectrum in terms of the effective field theory that is 
extended by additional fermionic species. Moreover, this result is completely independent 
of any random matrix theory representation of the quantum field theory. 

As noted in |]14|, the relationship ( |5.40| ) gives a much more compact form for the 



spectral density than that found in |T^. Here we have found an even more convenient 
expression for it, based on a coset parametrization of the finite volume gauge theory 
partition function. We stress once more that in the ordinary, unquenched case there is 
no need for this coset analysis because the group integral extends over U{Nf) up to a 
unitary group volume factor. However, in the supersymmetric case we are forced to deal 
directly with the coset space because this volume factor vanishes [^. The results of 
section complete the goal that was set out in section 1 of this paper, namely an analytical 
derivation of the QCD3 Dirac operator spectrum directly from quantum field theory. 
As a bonus, the supersymmetric form of the finite volume partially quenched partition 
function has thereby provided a new and much simpler expression for the microscopic 
spectral density. This illustrates the power behind the supersymmetric method. 
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6 Higher Order Spectral Correlation Functions 

Just as the equivalence between the universal random matrix theory and low energy 
effective field theory partition functions is not sufficient to establish the computability 
of the microscopic Dirac operator spectrum in random matrix theory, neither is merely 
the computation of the spectral one-point function. To complete the proof, one needs to 
extend the calculations to derive the generic /c-point spectral functions ps{ui, . . . , Uk, {^i})- 
It is clear that in this case one needs to consider a partially quenched quantum field theory 
involving N^ > k species of valence quarks, in which case the higher fc-point spectral 
correlators may be computed as the discontinuity across the cut of a /c-th order fermion 
susceptibility as 



27riSo/ ^-o+;Vi .,-^1 



(6.1) 



where 



E.(,^,.....^.i{.,}) = (-A) __^ 
In this section we will point out that the method described in the previous section can 



X 



• (6.2) 



be used to derive the function (|0|) from a partially quenched field theory partition func- 
tion involving A^^ = 2k species of valence quarks and supersymmetric coset G{nf;k,k) 
in ( |3.26| ). Thus in three spacetime dimensions, the program for computing the full mi- 
croscopic spectrum of ip may be completed in a straightforward way, unlike the four 
dimensional case. 

The calculation proceeds in exactly the same manner as in the previous section and 
we will therefore be very brief, only highlighting the salient points. The spectral fc-point 
function (|6.1| ) may be alternatively derived from |^ 



"^^ d f Ms 



k 



jLV/^r " viVE 






(6.3) 



The partially quenched field theory partition function in ( |6.3D may be evaluated in exactly 
the same way as described in section 5. For the regular part we find 

ZNf,2k{-^'i {f^2l-l, —fJ'2hP'2l-l, —^21} 
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det 



=1 1=1 ("^1 - l^2i-i){rni - /i2i-i) i<^«'<fc 



A^So(Ai2Z-l +/i2/') 



X det 



sin NTjq {rrii + rrij) sin A^Sq (mj + fi2i' 



sin A^So (yU2/- 



771 ,• 



A^So(mi + /i2i') 



(6.4) 



iVSo (yU2«-l + rrij) iVSo (yU2/_l + yU2/') 

where the second determinant is of dimension {nj + k) x [nj + k), while for the Efetov- 



Wegner term which yields the normalization (O) we find 

^^^2fc(-^; {f^2l-U -f^2hfi2l-l, -fi2l}) 
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, [sin NTiQ (nii + rrij] 

C[Qf^ - — 

^<i,j<nf NTiQ {iTii + TJlj) 



(6.5) 



The microscopic spectral fc-point function now follows from differentiating (|6.4| ) with 
respect to the valence quark masses as prescribed by (|6.3|) . Note that if we expand the 
kx k determinant in ( |6.4| ) as a sum over elements of the symmetric group 5*2^, then only 
the contribution from the identity permutation survives in the limit Jx2i-i + ^21 -^ 0, since 
it is only that term which cancels the factor of e'' in ( |6.3|) . In this way we arrive at 



det 



Psiui, ■■.,Uk] {uJi}j 



smh.{uji + ujj) sinh(ci;j + iui' 



uJi + ujj ijJi + iuv 

sinh(u;j — iui) mi{%ii — uy 



ijjj — lUi 



Ui - Ui' 



det 



sinh(u;j + ojj 

UJi + UJj 



(6.6) 



where the determinant in the numerator of (O) is of dimension [uf + /c) x (rij + k). For 
example, in the simplest case of quenched fermions, Nf = 0, the expression ( |6.6| ) reduces 
to 



psiui,...,Uk) = det 



1 sm{ui — Uy 

-K Ui- Uy 



(6.7) 



which is the expected representation in this case of the fc-point function in terms of the 
spectral sine-kernel of the unitary ensemble |]T1 . 

The expression 



is again an explicit representation of the formula |]TB| for the 
spectral fc-point function in terms of a ratio of two finite volume partition functions, one 
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of which now involves 2k additional species of fictitious quarks of imaginary mass, as 

k k nf 






u' 



1=1 i=l 



X — -^-^ 7 X ■ 6.8 



However, the representation (|6.6| ) is a new and much more explicit form of the microscopic 
fc-point spectral function. The equivalent representation of this function in terms oi akxk 
determinant of the random matrix theory spectral kernel |1^, |T^ leads to consistency 
conditions on the finite volume partition functions [0. The very explicit determinant 



formula ( |6.6| ) may provide a straightforward proof of this theorem, and hence of the 
equivalence of this field theoretical expression with that obtained from random matrix 
theory. We will not do so here, but we have checked this equivalence in a number of 
cases. As discussed in [0, these consistency conditions are related to the fact that the 
finite volume gauge theory partition function, which is an Itzykson-Zuber integral, is the 
r-function of an integrable KP hierarchy. 
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Appendix A Supersymmetric Limit 

Starting from the original random matrix theory partition function (|3.1| ) and following the 
exact same steps used in section 3.1 to arrive at ( |3.16|) , we may infer the supersymmetric 
representation 

ZN„N.{{m,};{fi,,fi,})= J DAe-^^T^^' SDet ^ (A - ^ >!') . (A.l) 

gl{Nf+N^\N^) 

In this appendix we will start by formally proving that the partition function ( |A.1| ) reduces 
to the expected one for ordinary QCD3 in the supersymmetric limit where /ij = /ij for 
each i = 1, . . . , A'^. For this, we decompose the {Nf + Ny\Nv) supermatrix A into the 
(1|1) supermatrix 

V Xi,Nf+i «Aii J 
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the 2- {Nf + 2N^ - 2) supervectors 



/ A 

'A = 



^. - I ^^^.^/+1 



^A,Nf+N^+l ^ 



*A = (A7v^.+i,a , ANf+N,+i,A) for Aj^Nf + l,Nf + N^ + l (A.3) 

of dimension (1|1), and the remaining {Nf + N^ — 1\N^ — 1) supermatrix which we denote 
by Ar. We may then write the integrand F(A) = e — 2^ ^^'^'' SDet^(A -iM') of {\AA\) 



as 



FiA) = F{A;{^A,^A};Ar) , (A.4) 

and in the degenerate case /Zj = fii it possesses the invariances 

f(uAU~'; {U^a, ^aU-'}; A,) = f(A; {^„, ^ J; A,) (A.5) 

for all U G G'L(lll). The partition function may then be written as 
^iv,,iv„({m,};{/i„^,}) = J DKH J D^aD^aZ{{'^b,^b};A,) 

(A.6) 
where 

z({^a,^a};A,)= j DAf(A;{^a,^a};A,) . (A.7) 

Using the invariance of the Haar measure one finds that ( [A.7| ) is an invariant function of 
the supervectors ( |A.3| ), 



z({U^A.^AU-'y,AA=z({^A,^A}\AA , UeGLil\l). (A.8 



It then follows from the Parisi- Sour las reduction ^7j that 

^7v,,iv„({m,};{/i„/i,}) = J DA,z[{0,Oy,A,). (A.9) 

gl(Nf+N^-l\N^-l) 

Next we may invoke the Efetov-Wegner theorem |2^ which states that 

/ DA G{A) = G(0) (A.IO) 

9«{1|1) 

for any supersymmetric invariant function G, i.e. G{A) = G{UAU^^) for all U G G'L(lll). 
Applying this result to the integral ^({0, 0}; Ar), we arrive at 

^iv„jv.({m,}; {/.,,,/.,}) = j DA, f(0; {0, 0}; A^) . (A.ll) 

gliNf+N^-l\N.u-l) 
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Now we repeat this procedure by decomposing Ar analogously, starting with a (1|1) 
supermatrix Aj. with diagonal bosonic elements A7v^+2,Ar^+2 and zA22, as in ([A.2| ). We 
iterate this reduction until all the Grassmann components and the Ny x N^ boson-boson 
block of the original supermatrix A have been eliminated. The final result is the reduced 
partition function 



ZNf,N^({'mi};{fii,fii}) = Zn^.^q{M) = I DX e 

uiNf) 



— ^'^ det 



N 



(^X-iM) . (A. 12) 



It is independent of the valence quark masses and is the standard one for ordinary QCD3 
expressed as an integral over the physical fiavour space [0. In the microscopic limit, 
it becomes the finite volume partition function ( |1.6| ), Zj^'^q{A4) = ZjfiM.) [^0. In 
the original field theory formulation, this reduction can be understood in perturbation 
theory by noting that each bosonic superpartner to the valence fermions contributes a 
Feynman diagram of equal magnitude but opposite sign to the partition function. This 
cancellation of graphs with valence fermion loops is precisely what is required for the 
associated spectral density to be equal to the QCD3 spectral density. Note that this 
normalization is non-trivial, since the invariances of the function G{A) would naively 
imply that the integral (|A.1CI|) vanishes because of the Grassmann integrations. Indeed, 
the right-hand side of (|A.1CI|) is an example of an Efetov-Wegner boundary term which is 
characteristic of superintegrals |2^-|2^ . 

A similar argument can be used to show that the supersymmetric form (|3.16|) of the 
random matrix theory partition function ( p.3|) reduces to (|A.1|) in the limit TW ^ 00. For 
this, we write the superdeterminant in ( |3.16| ) as 



SDet A-iM 



Nf . 

li- 

1 rrij 
j=i 3 



SDet 



dmg{\Nj+N.\M \liv„) A 



+ diag(A1,/il,...,/iAr^ lAr^,/il,...,/iAr„) . (A.13) 



We now use ( |3.17] ) to expand ( |A. 13| ) for ttIq, -^ 00 in the Nf x Nf boson-boson block 
matrix corresponding to the physical quark superpartners. In this limit, the argument of 
the superdeterminant in ( A.13| ) is a GL{Nf, C) invariant function of this block, and also 
of the 2A^„ supervectors in C ^' ^ which comprise the Nf x A^^ and N^ x Nf blocks of A 
and their supersymmetric counterparts. By applying the above reduction theorems, we 
thereby find that (|3.16|) reduces to ( |A.1| ) in the large mass limit. The question of whether 
the supersymmetric and large mass limits commute with the thermodynamic limit N ^ 00 
is somewhat more subtle. The appropriate large mass expansion for A^ — »• cx) is described 
in section 6, along with the pertinent Efetov-Wegner terms for the supersymmetric limit. 
For a more detailed proof that the limits A^ ^ 00 and rha -^ C)0 are commutable at 
least for the purpose of computing spectral correlation functions, see [^, where a similar 
supersymmetrization using fictitious supersymmetric sea quark partners was applied to 
the partition function of the chiral Gaussian unitary ensemble relevant for QCD4. 
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Appendix B Supersymmetric Itzykson-Zuber For- 
mula 

In this appendix we will give a simple derivation of the supersymmetric extension of the 
Itzykson-Zuber formula for the unitary supergroup U{N\M) |2^. The integral is 



X[X,Y-k]= I z^f/ e'^STt(^^^^0 (B.l) 

U{N\M) 

where 

N M 

DU = 11 dU,, dU*^ 6 (Ea U^A U;^ - 5,,) n dUap dUlp 5 (Ea UaA U;^ - 6^p) 

i,j=l a, 13=1 

N M r. o N+M Af A/ Q Q N+M 

is the invariant Haar-Berezin measure on U{N\M), k G C is a constant parameter, and 
X, y G u{N\M). The super-unitary invariance of the Haar measure implies that X 
and Y may be taken to be diagonal without any loss of generality. Let xa = {xi,Xa) 
and i/A = {yi,ya) be their respective supereigenvalues, where here and in the following 
capital Latin letters A = 1, . . . ,N + M label the complete set of supersymmetric indices, 
while lower case Latin letters i = 1, . . . ,N label the fermionic indices and Greek letters 
a = 1,...,M the bosonic indices (This is the same notation used in the text). We 
denote by e{A) the Grassmann grading of the index A defined by e{i) = 1 mod 2 and 
e{a) = mod 2. 

The crucial observation that we shall make here is that the integral ( |B.1| ) is effectively 
defined over the homogeneous superspace U{N\M)/U{1)^^^^ . In the purely bosonic case 
M = 0, such an integral would be over a coadjoint orbit of the ordinary unitary group and 
it would define a dynamical system which satisfies the hypotheses of the Duistermaat- 
Heckman theorem (see |3^ for a detailed exposition of the subject). This means that 



the semi-classical approximation to the integral, obtained by summing over all extrema 
(minima, maxima and saddle-points) of the argument of the exponential, is exact. In the 
present case, we need the analog of the Duistermaat-Heckman theorem for supermani- 
folds. The conditions under which the stationary phase approximation is exact for such 
superintegrals are discussed in |3^. In the following we will assume that these criteria 



are met by the integral ( |B.1| ) and simply evaluate it in the saddle-point approximation. 
This is justified by the fact that the ordinary integration manifolds supporting the unitary 
supergroup are symmetric spaces, in the usual sense. The basic point is that the superin- 
tegral (p.l| ) contains the same, large amount of symmetries that its bosonic counterpart 
has, which is the feature that is always responsible for the exactness of the semi-classical 
approximation in such instances [^ . 
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For this, we first need to find the extrema of the function 

N+M 



A,B=1 



n[U] = STr (Xf/Ff/t) = Yl {-ly^^^'-'xAyB Uab 
By using the identity 



d 



dU. 



AB 



■uln = i-if'''-'^'''^'ul^u 



BD 



for U G U{N\M), we find that the saddle-point equation reads 



X , UYU^ 







which for X and Y diagonal becomes 

N+M 

XA Y. UacVcUIb = (-1^ 

C=l 



N+M 



xb J2 UacVcUIb 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



c=i 



Up to irrelevant elements of the Cartan subgroup U{1)^^^'' of U{N\M)^ the only so- 
lutions U of (p.6| ) are those matrices which permute the eigenvalues of the matrix F, 
i.e. YjcUacVcUcb = yp{A)5AB with P G Sn+m- It is easy to see that there are no 
Grassmann-odd permutation matrices P G U{N\M) that can map bosonic and fermionic 
indices into one another, i.e. for which i = P[a). Therefore, only Grassmann-even P can 
occur and they take the generic form 



n 







AB 



n„;3 



with u e Siv , n G Sm 



>M ■ 



(B.7) 



The saddle-point approximation thereby dictates to sum over all elements of the discrete 
Weyl subgroup of the ordinary Lie group supporting U{N\M). 

We now set f/ = (11 © 11) e*^ in (p.3| ), with L G u{N\M) an infinitesimal Hermitian 
supermatrix, and expand the function 7i[f/] to quadratic order in L. The semi-classical 
approximation to the unitary supermatrix integral ( [B.l|) , obtained by summing over all 
extrema, thereby reads 

1 / N M \ 

x[x,Y;k] = — — Y. J2 ^^P^\J2^iym)-Y.^c.yn{a)] 



TieSN ueSM 

N 



X 



d d 



Af „ M N M 

/ n dL.^ / n dL^p ® n n .. .. 

■j) {yui^) - yu{j)) 

1 iw 2 

- 22 |-^a/3 \^a — Xf3J {yn{a) " l/n(/3)^ 



u{N) 

X exp K 



u{M) 
1 ^ I 



lAj li lAJ -l 



o,/3=l 

N M 



+ H II I ^»" (^« ~ ^o) (yn(i) - yii{a) 



i=l a=l 



(B.8) 
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By evaluating the Gaussian integrals over the complex bosonic variables Lij, i ^ j and 
Lapi Q- ^ I3i the real bosonic variables La and Lqq,, and the complex Grassmann variables 



LjQ,, we arrive at 



X[X,Y-k] 



1 



N 



N\M\ 



E E n 



,i<-xtyii{i) 



nG5jv ne5A./ »=i 
sgnll 



N(N-l) 

"" \~) A[a:]A[y] 



M 

n 



8q yn(a) 



M(M-l) 
2 



sgnll 

A[i]A|j 



K 



ATM 



A[x, x]A[y,y] 



where 



A [A] 



det 



Ar 






(B.9) 



:b.io) 



is the Vandermonde determinant, and 



A 



A, A 



N M 

nn(A,-A. 

i=l a=l 



(B.ll) 



In arriving at ( |B.9| ) we have used the properties 

n (^n(.) - Ano)) = sgn(n)A[A] 



N M 



nn(An 



(i) 



\n(a) 



i=l Q=l 

Summing over the permutations in 



A 



A, A 



(B.12) 



then leads to 



I[X,Y;k] 



N(N-l) M(M-l) 
(ZTT) 2 "^ 2 JV(iV-l) 
K 2 



Af(M-l) 



iV!M! 



X A[a;,a;] A[?/,y] 



-K 



det 



CjJ/j 



det 

o,/3 



-KXaVfs 



A[x]A[ 



A[x]A[y] 



(B.13) 



which is the standard supersymmetric generalization of the Itzykson-Zuber formula [22 



The above derivation can be extended in a straightforward fashion to the generaliza- 
tion of (|B.1|) over any connected, compact, semi-simple Lie supergroup G, with X and 
Y elements of the Cartan subalgebra of the Lie superalgebra of G. By extending the 



formalism described at length in [^, one may in this way derive the appropriate super- 
symmetric extension of the Harish-Chandra formula [^ in the form conjectured in p3| . 
The present method gives a much simpler and compact way of deriving these supersym- 
metric integration formulas, in contrast to the superspace heat kernel and supergroup 
theoretic methods employed in [^ . 



^The expression ( |B.13| ) agrees with those obtained in |g^ up to the overall numerical prefactor which 
is related to the volume of the ordinary unitary group supporting U{N\M) in the bosonic Haar measure 
and also the volume of the corresponding Weyl subgroup. 
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